Test of Fermi gas model and plane-wave impulse approximation 
against electron-nucleus scattering data 
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A widely used relativistic Fermi gas model and plane-wave impulse approximation approach are 
tested against electron-nucleus scattering data. Inclusive quasi-elastic cross section are calculated 
and compared with high- precision data for 12 C, 16 0, and 40 Ca. A dependence of agreement between 
calculated cross section and data on a momentum transfer is shown. Results for the 12 C(u^, fi~) 
reaction are presented and compared with experimental data of the LSND collaboration. 
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I. INTRODUCTION 

A precise description of neutrino-nucleus (vA) interac- 
tions in intermediate energy region is important for the 
interpretation of measurements by many neutrino exper- 
iments. The understanding of their sensitivity to neu- 
trino properties, evaluation of the neutrino fluxes and 
spectra depend on the accuracy to which the vA cross 
sections are known. This is in particular crucial in analy- 
sis of the long-base line neutrino oscillation experiments 
in which the parameter of neutrino oscillation Am 2 is 
determined using the total number of detected events 
and the distortions in the energy distribution of the de- 
tected muons caused by neutrino oscillation. On the 
other hand the neutrino-nucleus cross sections contain 
contributions from both axial-vector and vector currents 
and thus provide complementary information to that pro- 
vided by electron-nucleus scattering, which is sensitive 
only to the nuclear vector current. 

In many experiments the neutrino fluxes in sub-GeV 
and GeV energy region are used. At such energies the 
charged-current quasi-elastic (QE) neutrino-nucleus in- 
teractions give the main contribution to the detected 
events. Sizable nuclear effects have been observed in lep- 
ton scattering off nucleus at energies less than a few GeV. 
They indicate that the nuclear environment plays an im- 
portant role even at energies and momenta larger than 
those involved in typical nuclear ground state processes. 
The understanding of nuclear effects is relevant for the 
long-base line neutrino experiments in order to control 
the corresponding systematic uncertainties. 

Many Monte-Carlo (MC) [l| codes developed for sim- 
ulation of the neutrino detectors response are based on 
a simple picture of a nucleus as a system of quasi-free 
nucleons, i.e. a relativistic Fermi gas model (RFGM) 
0. It takes into account Fermi motion of nucleons in- 
side the nucleus and Pauli blocking effect. Unfortunately 



the uncertainties of a few vA scattering data at low and 
intermediate energies don't allow us to estimate the ac- 
curacy of this model. On the other hand, as follows from 
vast high-precision electron scattering data the RFGM 
neglects some important nuclear effects. So, the model 
which is used for calculation of neutrino scattering off 
nucleus should first be tested against electron scattering 
data. 

In the present work we intent to evaluate the accu- 
racy of the Fermi gas model and plane-wave impulse 
approximation (PWIA) [3,4,5]. The electron QE cross 
sections are calculated in the framework of these mod- 
els and compared with high-precision data for different 
nuclei. This comparison shows that the agreement be- 
tween predictions of these models and data depends sig- 
nificantly on the momentum transfer to the target. We 
applied the RFGM and plane-wave impulse approxima- 
tion to 12 C(f M ,^~) reaction also. 

The formalism of an inclusive charged current lepton- 
nucleus QE scattering is given in Sec. 2. Results are pre- 
sented and discussed in Sec. 3 and some conclusions are 
drawn in Sec. 4. 



II. FORMALISM OF THE INCLUSIVE 
QUASI-ELASTIC SCATTERING 

In electromagnetic and weak charge current process 
electrons (neutrinos) interact with nuclei via the ex- 
change of photons or W-boson and charged leptons are 
produced in the final state. In an inclusive reaction, in 
which incident electron (a el ) or neutrino (cr cc ) with four- 
momentum ki — (si, kj) is absorbed by nucleus with mass 
rriA and only the out-going lepton with four-momentum 
kf = (e/,k/) and mass mi is detected, the cross section 
is given by contracting the lepton tensor and the nuclear 
tensor 
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where a ~ 1/137 is the fine-structure constant, G ~ 
1.16639 xlO" 11 McV~ 2 is the Fermi constant, 6 C is the 
Cabbibo angle (cos6» c «0.9749), Q 2 = -q 2 = (fc; - k f ) 2 , 
and q = (w, q) is the four-momentum transfer. 

The lepton tensor can be written, by separating the 
symmetrical Ig 1 ' and antisymmetric components V% as 
follows 
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where M is the nucleon mass and u^ v — i[y^ v }/2. Fy 
and Fm are the isovector Dirac and Pauli nucleon form 
factors, taken from Ref.[6]. Fa and Fp are axial and 
pseudo-scalar form factors, parametrized as 
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where Fa{0) — 1.267, is pion mass, and Ma — 1.032 
GeV is axial mass. 

The general covariant form of the nuclear tensors is 
obtained in terms of two four-vectors, namely the four- 
momenta of target p M and q 11 . The electromagnetic and 
charged-current nuclear tensors can be written as 



= 2{k<tk) + tfk? - gTh ■ k f ), 
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where e^ va ^ is the antisymmetric tensor with e 0123 = 
— eoi23 = I- Assuming the reference frame in which the 
z-axis is parallel to the momentum transfer q = — ky 
and the y-axis is parallel to k^ x k/, the symmetrical 
components l s 3 : ,lg l ' ,l Z g and the anti-symmetrical ones 
l A x ,l x A z ,l A z , as well as those obtained from them by ex- 
changing their indices, vanish. 

(el) 

The electromagnetic W^u and weak charged-current 

Wpif* hadronic tensors are given by bilinear products of 
the transition matrix elements of the nuclear electromag- 



netic (weak charged current) operator J ; 
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between 



the initial nucleus state \A) of energy E and final state 
\Bf) of energy Ef as 

W$ ){cc) = ^{B f \J { f ){cc) \A) x (A\J^ cc ^\B f ) 
f 

x5{E +u)-E f ), (6) 

where the sum is taken over the undetected states. 

The transition matrix elements are calculated in the 
first order perturbation theory and in impulse approxi- 
mation, i.e. assuming that the incident lepton interacts 
with the single nucleon while other ones behave as spec- 
tators. The nuclear current operator jjf l ^ cc \q) is taken 
as the sum of single- nucleon currents (q), i.e. 
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where Wi are nuclear structure functions which depend 
on two scalars Q 2 and p ■ q. Form Eqs.(10),(ll) it fol- 
lows that in our reference frame W 0x , W° v , and W xz , 
vanish together with tensor components obtained from 
them by exchanging their indices. Therefore, obtained 
from contraction between lepton Eqs.(3),(4),(5) and nu- 
clear Eqs.(10),(ll) tensors, the inclusive cross sections for 
the QE electron (neutrino)-nucleus scattering are given 
by 



da^ 



defdfl 



— Pm{V l ri L + v T K T ) 



(el) r>( e 'h 



(12) 



da (cc ^ 
defdfl 



G 2 cos 2 0, 



kfEf^oR^ + vpRp 



? (cc) 



(27T) 2 

-volR^ + vllR^ ± v^Ri^), (13) 



a 2 cos 2 (fl/2) 
4s 2 sm A (6/2)' 



(14) 



where <jm is the Mott cross section and 9 is the lepton 
scattering angle. The coefficients and v are obtained 
from lepton tensors and can be written as 
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FIG. 1: Nucleon momentum distribution corresponding to Eq.(28) (solid lines) and Eq.(26) (dotted lines). The momentum 
distribution no is given by short-dashed line. The open squares represent results obtained in Ref. [8] . The full triangles represent 
the values of no(p) obtained in Ref. [9]. 
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where k = |k/|/e/. The nuclear response functions 
and i?( cc ) are given in terms of components of nuclear 
tensors as follows [4], [7] 
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In order to evaluate nuclear response functions we con- 
sider the RFGM and PWIA approach based on assump- 
tion that the virtual photon interacts with off-shell nu- 
cleon and neglecting interaction of the knocked out nu- 
cleon with the residual nucleus. In the PWIA the four- 
momenta of the initial nucleus A, the bound off-shell nu- 
cleon N, and the final state B are 
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respectively. Here m* B 
mass and intrinsic energy of the final (A-l)-nucleon state, 
respectively. Within the above assumption the nuclear 
structure functions can be written as follows 
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+similar terms for the neutrons, 
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Here, Z is the number of protons, W"'° Tt (N — p, n) are 



foff {N = , 

the off-shall nucleon structure functions that are given 
in the terms of nucleon form-factors. S N (\p\,E) is the 
nucleon spectral function and kinematic factors Cjj , Dij 
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FIG. 2: Comparison of theoretical and experimental cross sections for 12 C. The data are taken from Refs.[10] (filled circles), 
[11] (filled squares), [12] (filled triangles), [13] (open circles), [14] (open squares), and [15] (stars). 
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FIG. 3: Comparison of theoretical and experimental cross sections for 1 O. The data are taken from Refs.[ll] (filled circles), 
and [16] (filled triangles) . 
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where cost = p • q/|p • q|. In this paper we assume that 
W^'° are identical to free nucleon structure function 
. The parametrization of W? 1 is taken from Refs. [2,3] 
as follows 
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where Pq = [(p + q) 2 + M 2 ] 1 / 2 is energy of the knocked- 
out nucleon. 

The nucleon spectral function S N (\p\,E) in the PWIA 
represents probability to find the nucleon with the mo- 
mentum p and the removal energy E in the ground state 
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FIG. 4: Comparison of theoretical and experimental cross sections for 40 Ca. The data are taken from Ref.[18]. 



of the nucleus. In the commonly used Fermi gas model, 
that was described by Smith and Moniz [2], nucleons in 
nuclei are assumed to occupy plane wave states in uni- 
form potential while the knocked-out nuclcon is outside 
of the Fermi sea. The Fermi gas model provides the sim- 
plest form of the spectral function which is given by 

s fg (e,\ p \) = -J^e^-ipDeflp + qi-p,.) 

xS[(\p\ 2 + M 2 ) 1 / 2 - e b - E)], (27) 
where pf denotes the Fermi momentum and a parame- 



ter e& is effective binding energy, introduced to account 
of nuclear binding. The QE lepton-nucleus reactions are 
complicated processes, involving nuclear many body ef- 
fects. The calculation of the nuclear spectral function for 
complex nuclei requires to solve many body problem. In 
this paper we consider also a phenomenological model us- 
ing PWIA approach with the spectral function which in- 
corporates both the single particle nature of the nucleon 
spectrum at low energy and high-energy and high mo- 
mentum components due to NN-correlations in ground 
state. Following [4,5] we separate the full spectral func- 
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FIG. 5: Comparison of theoretical and experimental cross sections for Ca. The data are taken from Refs.[10] (filled circles), 
and [18] (filled triangles) . 



tion into two parts 

S(E,p) = S (E,p) + S 1 (E,p). (28) 

The integration of Eq. (27) over energy gives nucleon mo- 
mentum distribution, 

n(p) = J ^S(E,p) = no(p)+ni(p). (29) 

The spectral function is normalized according to 

J d -^S(E, P) = 1. (30) 

The low-energy part (So) can be approximated by the 
following expression 

5 (£,p) = 2™ (pM(£-£ (1) -p 2 /2M A _ 1 ), (31) 

with the nucleon separation energy averaged over 
residual configurations of A — 1 nucleons with low excita- 
tion energies and recoil energy p 2 /2Ma_i, no(|p|) is the 
corresponding part of nucleon momentum distribution. 

The high-excitation part (Si ) of the spectral function 
is determined by excited states with one or more nucle- 
ons in continuum. The detailed description of this model 



is given in Rcfs.[4,5] as well as parametrization of rio(p) 
and ni(p), which fit the result of many-body calcula- 
tions of nuclear momentum distribution. As follows from 
these calculations the low momentum part incorporates 
about 80% of the total normalization of spectral func- 
tion, while the other 20% are taken by the high momen- 
tum part. The nucleon momentum distributions n(|p|) 
and ufg(\p\) are shown in Fig.l. The normalization of 
n(p) and Ufg(p) is / dpp 2 n(p) = 1, where p = |p|. The 
distributions Ufg for various nucleus 12 C, 16 and 40 Ca 
were calculated using the value of parameters pf = 221 
MeV, e b = 25 MeV ( 12 C); p F = 225 MeV, e b = 27 MeV 
( 16 0); and p F = 249 MeV, e b = 33 MeV) ( 40 Ca) [10]. 

III. RESULTS 

There is vast high-precision data for electron scatter- 
ing off nucleus 12 C, 16 0, and 40 Ca. Hence these nuclei 
are taken at the focus of the present work. Data on in- 
clusive cross sections for a number of nuclei (A between 
6 and 208) but only one set of kinematics were obtained 
early in Ref. [10]. Carbon data are available from experi- 
ments [11]-[15]. For oxygen target the experiments were 
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FIG. 6: Comparison of theoretical and experimental cross sections for Ca. The data are taken from Refs.[17] (open squares), 
and [18] (filled triangles). 



performed by SLAC [11] and Frascaty [16] groups. For 
calcium target the inclusive cross section have been mea- 
sured in experiments [10], [17]-[19]. All these data were 
used in our analysis. 

Using both the relativistic Fermi gas model and the 
PWIA approach described above, we calculated the in- 
clusive cross sections for given kinematics (energies and 
angles) and compared them with data. The results are 
presented in Figs. 2-6 for 12 C, 16 0, and 40 Ca respectively. 
The solid lines are the results in the Fermi gas model, 
while short-dashed lines are results in the PWIA. The dif- 
ferences can be seen from these figures in which the cross 
sections as functions of u) or invariant mass produced 
on a free nucleon W are plotted. At the maximum of 
the cross sections both models overestimate the measured 
values. We evaluated the differences between predicted 
(a ca i) and measured {a data) quantities A = o calc - a data . 
A(|q|) as a function of three-momentum transfer |q|, is 
shown in Fig. 7, from which it is clear that the A(|q|) 
decreases with |q| from 30-^50% at |q| < 200 MeV to 
10-=-15% at |q| > 500 MeV. 

In Refs. 01 1 transverse and longitudinal 

response functions have been extracted for 200 



MeV< |q| < 500 MeV . It has been shown that the 
RFGM overestimates the observed longitudinal response 
for about 40% [13] (~20% 18]). At low |q| this model 
also overestimates the magnitude of the transverse re- 
sponse function. At higher values of |q| it do better at 
reproducing the value of the ■ 

The predictions of both models are compared with 
the experimental result of the LSND collaboration at 
Los Alamos for 12 C(z^,/i~) reaction |2(|. The calcu- 
lations are flux-averaged over the Los Alamos neutrino 
flux. The mean energy of neutrino flux above thresh- 
old is 156 MeV. The comparison is shown in Fig. 8 where 
the calculated muon energy distributions are normalized 
to the experimental total number events. We note that 
both models do not give an accurate descriptions of the 
shape of the muon spectrum. The flux-averaged cross 
section integrated over the muon energy is 17.8xl0~ 40 
cm 2 in the case of the RFGM and 26.8xl0" 40 cm 2 in the 
PWIA. The experimental value is (10.6±0.3±1.8)xl0^ 40 
cm 2 . The result obtained by other calculation in the 
framework of the Fermi gas model with local density 
approximation [21] gives also larger value er=(16.7 ± 
1.37) xlO~ 40 cm 2 . 
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function A(|q|) decreases with three-momentum transfer 
from 30^-50 % at |q| < 200 MeV to 10^-15 % at |q| > 
500 MeV. Therefore these models overestimate also the 
cross sections at low Q 2 = |q| 2 — u> 2 . 

We applied the RFGM and PWIA approach to 
12 C(i/ /i , /i~) reaction. The flux-averaged total cross sec- 

300 




FIG. 8: The distribution of muon kinetic energy for inclusive 
12 C (Vfj,,fj,~) reaction. Experimental data from Ref.[20]. The 
results of the RFGM (solid linen histogram) and the PWIA 
approach (short-dashed line histogram) are normalized to the 
data. 



FIG. 7: Differences between calculated and measured values 
of cross sections at maximum for 12 C, ie O, and 40 Ca as func- 
tions of three-momentum transfer |q|. The filled triangles 
correspond to the Fermi gas model results and open circles 
correspond to the PWIA approach. 

IV. CONCLUSIONS 

In this work we have tested the widely used relativis- 
tic Fermi gas model and plane-wave impulse approxima- 
tion approach against electron-nucleus scattering data. 
We calculated the inclusive QE cross sections and com- 
pared them with high-precision data for 12 C, 16 O, and 
40 Ca in a wide region of incident energy and momentum. 
We evaluated the differences A between predicted and 
measured QE cross section at the maximum and found 
that both models overestimate the measured values. The 



tions and muon energy distributions were calculated and 
compared with experimental results of the LSND collab- 
oration. The calculated cross sections are significantly 
larger than the experimental ones and both models do 
not give an accurate description of the shape of muon 
spectrum. 

In conclusion we note that the inclusion of final state 
interaction effects along with realistic spectral function 
may significantly correct the description of the data at 
low momentum transfer, as was pointed in Ref. [16]. 
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